Spatially resolved single-cell atlas of the lung in fatal Covid19 in an African population reveals a distinct cellular signature and an interferon gamma dominated response (preprint)

Postmortem single-cell studies have transformed understanding of lower respiratory tract diseases (LRTD) including Covid19 but there is almost no data from African settings where HIV, malaria and other environmental exposures may affect disease pathobiology and treatment targets. We used histology and high-dimensional imaging to characterise fatal lung disease in Malawian adults with (n=9) and without (n=7) Covid19, and generated single-cell transcriptomics data from lung, blood and nasal cells. Data integration with other cohorts showed a conserved Covid19 histopathological signature, driven by contrasting immune and inflammatory mechanisms: in the Malawi cohort, by response to interferon-gamma (IFN-{gamma}) in lung-resident alveolar macrophages, in USA, European and Asian cohorts by type I/III interferon responses, particularly in blood-derived monocytes. HIV status had minimal impact on histology or immunopathology. Our study provides data resources and highlights the importance of studying the cellular mechanisms of disease in underrepresented populations, indicating shared and distinct targets for treatment.

The SARS-CoV-2 Omicron sub-variant BA.2.86 is attenuated in hamsters (preprint)

SARS-CoV-2 variants have emerged throughout the COVID-19 pandemic. There is a need to risk-assess newly emerged variants in near "real-time" to estimate their potential threat to public health. The recently emerged Omicron sub-variant BA.2.86 raised concerns as it carries a high number of mutations compared to its predecessors. Here, we assessed the virulence of BA.2.86 in hamsters. We compared the pathogenesis of BA.2.86 and BA.2.75, as the latter is one of the most virulent Omicron sub-variants in this animal model. Using digital pathology pipelines, we quantified the extent of pulmonary lesions measuring T cell and macrophage infiltrates, in addition to alveolar epithelial hyperplasia. We also assessed body weight loss, clinical symptoms, virus load in oropharyngeal swabs, and virus replication in the respiratory tract. Our data show that BA.2.86 displays an attenuated phenotype in hamsters, suggesting that it poses no greater risk to public health than its parental Omicron sub-variants. Article summary lineThe newly emerged Omicron sub-variant BA.2.86 is attenuated in hamsters.

Disease trajectories in hospitalized COVID-19 patients are predicted by clinical and peripheral blood signatures representing distinct lung pathologies. (preprint)

Linking clinical biomarkers and lung pathology still is necessary to understand COVID-19 pathogenesis and the basis of progression to lethal outcomes. Resolving these knowledge gaps enables optimal treatment approaches of severe COVID-19. We present an integrated analysis of longitudinal clinical parameters, blood biomarkers and lung pathology in COVID-19 patients from the Brazilian Amazon. We identified core signatures differentiating severe recovered patients and fatal cases with distinct disease trajectories. Progression to early death was characterized by rapid and intense endothelial and myeloid activation, presence of thrombi, mostly driven by SARS-CoV-2 + macrophages. Progression to late death was associated with systemic cytotoxicity, interferon and Th17 signatures and fibrosis, apoptosis, and abundant SARS-CoV-2 + epithelial cells in the lung. Progression to recovery was associated with pro-lymphogenic and Th2-mediated responses. Integration of antemortem clinical and blood biomarkers with post-mortem lung-specific signatures defined predictors of disease progression, identifying potential targets for more precise and effective treatments.

Phenotyping the virulence of SARS-CoV-2 variants in hamsters by digital pathology and machine learning (preprint)

SARS-CoV-2 has continued to evolve throughout the COVID-19 pandemic, giving rise to multiple variants of concern (VOCs) with different biological properties. As the pandemic progresses, it will be essential to test in near real time the potential of any new emerging variant to cause severe disease. BA.1 (Omicron) was shown to be attenuated compared to the previous VOCs like Delta, but it is possible that newly emerging variants may regain a virulent phenotype. Hamsters have been proven to be an exceedingly good model for SARS-CoV-2 pathogenesis. Here, we aimed to develop robust quantitative pipelines to assess the virulence of SARS-CoV-2 variants in hamsters. We used various approaches including RNAseq, RNA in situ hybridization, immunohistochemistry, and digital pathology, including software assisted whole section imaging and downstream automatic analyses enhanced by machine learning, to develop methods to assess and quantify virus-induced pulmonary lesions in an unbiased manner. Initially, we used Delta and Omicron to develop our experimental pipelines. We then assessed the virulence of recent Omicron sub-lineages including BA.5, XBB, BQ.1.18, BA.2 and BA.2.75. We show that in experimentally infected hamsters, accurate quantification of alveolar epithelial hyperplasia and macrophage infiltrates represent robust markers for assessing the extent of virus-induced pulmonary pathology, and hence virus virulence. In addition, using these pipelines, we could reveal how some Omicron sub-lineages (e.g., BA.2.75) have regained virulence compared to the original BA.1. Finally, to maximise the utility of the digital pathology pipelines reported in our study, we developed an online repository containing representative whole organ histopathology sections that can be visualised at variable magnifications (https://covid-atlas.cvr.gla.ac.uk). Overall, this pipeline can provide unbiased and invaluable data for rapidly assessing newly emerging variants and their potential to cause severe disease.

Global landscape of the host response to SARS-CoV-2 variants reveals viral evolutionary trajectories (preprint)

A series of SARS-CoV-2 variants of concern (VOCs) have evolved in humans during the COVID-19 pandemic: Alpha, Beta, Gamma, Delta, and Omicron. Here, we used global proteomic and genomic analyses during infection to understand the molecular responses driving VOC evolution. We discovered VOC-specific differences in viral RNA and protein expression levels, including for N, Orf6, and Orf9b, and pinpointed several viral mutations responsible. An analysis of the host response to VOC infection and comprehensive interrogation of altered virus-host protein-protein interactions revealed conserved and divergent regulation of biological pathways. For example, regulation of host translation was highly conserved, consistent with suppression of VOC replication in mice using the translation inhibitor plitidepsin. Conversely, modulation of the host inflammatory response was most divergent, where we found Alpha and Beta, but not Omicron BA.1, antagonized interferon stimulated genes (ISGs), a phenotype that correlated with differing levels of Orf6. Additionally, Delta more strongly upregulated proinflammatory genes compared to other VOCs. Systematic comparison of Omicron subvariants revealed BA.5 to have evolved enhanced ISG and proinflammatory gene suppression that similarly correlated with Orf6 expression, effects not seen in BA.4 due to a mutation that disrupts the Orf6-nuclear pore interaction. Our findings describe how VOCs have evolved to fine-tune viral protein expression and protein-protein interactions to evade both innate and adaptive immune responses, offering a likely explanation for increased transmission in humans.

Enhanced innate immune suppression by SARS-CoV-2 Omicron subvariants BA.4 and BA.5 (preprint)

SARS-CoV-2 adaptation to its human host is evidenced by the emergence of new viral lineages with distinct genotypic and phenotypic characteristics, termed variants of concern (VOCs). Particular VOCs have become sequentially dominant globally (Alpha, Delta, Omicron) with each evolving independently from the ancestral Wuhan strain. Omicron is notable for its large number of Spike mutations found to promote immune escape and re-infection. Most recently, Omicron BA.4 and BA.5 subvariants have emerged with increasing levels of adaptive immune escape threatening vaccine effectiveness and increasing hospitalisations. Here, we demonstrate that the most recent Omicron variants have enhanced capacity to antagonise or evade human innate immune defenses. We find Omicron BA.4 and BA.5 replication is associated with reduced activation of epithelial innate immune responses versus earlier BA.1 and BA.2 subvariants. We also find enhanced expression of innate immune antagonist proteins Orf6 and N, similar to Alpha, suggesting common pathways of human adaptation and linking VOC dominance to improved innate immune evasion. We conclude that Omicron BA.4 and BA.5 have combined evolution of antibody escape with enhanced antagonism of human innate immunity to improve transmission and possibly reduce immune protection from severe disease.

The hyper-transmissible SARS-CoV-2 Omicron variant exhibits significant antigenic change, vaccine escape and a switch in cell entry mechanism (preprint)

Vaccination-based exposure to spike protein derived from early SARS-CoV-2 sequences is the key public health strategy against COVID-19. Successive waves of SARS-CoV-2 infections have been characterised by the evolution of highly mutated variants that are more transmissible and that partially evade the adaptive immune response. Omicron is the fifth of these Variants of Concern (VOCs) and is characterised by a step change in transmission capability, suggesting significant antigenic and biological change. It is characterised by 45 amino acid substitutions, including 30 changes in the spike protein relative to one of the earliest sequences, Wuhan-Hu-1, of which 15 occur in the receptor-binding domain, an area strongly associated with humoral immune evasion. In this study, we demonstrate both markedly decreased neutralisation in serology assays and real-world vaccine effectiveness in recipients of two doses of vaccine, with efficacy partially recovered by a third mRNA booster dose. We also show that immunity from natural infection (without vaccination) is more protective than two doses of vaccine but inferior to three doses. Finally, we demonstrate fundamental changes in the Omicron entry process in vitro, towards TMPRSS2-independent fusion, representing a major shift in the replication properties of SARS-CoV-2. Overall, these findings underlie rapid global transmission and may alter the clinical severity of disease associated with the Omicron variant.

The P681H mutation in the Spike glycoprotein confers Type I interferon resistance in the SARS-CoV-2 alpha (B.1.1.7) variant (preprint)

Variants of concern (VOCs) of severe acute respiratory syndrome coronavirus type-2 (SARS-CoV-2) threaten the global response to the COVID-19 pandemic. The alpha (B.1.1.7) variant appeared in the UK became dominant in Europe and North America in early 2021. The Spike glycoprotein of alpha has acquired a number mutations including the P681H mutation in the polybasic cleavage site that has been suggested to enhance Spike cleavage. Here, we show that the alpha Spike protein confers a level of resistance to the effects of interferon-{beta} (IFN{beta}) in lung epithelial cells. This correlates with resistance to restriction mediated by interferon-induced transmembrane protein-2 (IFITM2) and a pronounced infection enhancement by IFITM3. Furthermore, the P681H mutation is necessary for comparative resistance to IFN{beta} in a molecularly cloned SARS-CoV-2 encoding alpha Spike. Overall, we suggest that in addition to adaptive immune escape, mutations associated with VOCs also confer replication advantage through adaptation to resist innate immunity.

Mutations that adapt SARS-CoV-2 to mustelid hosts do not increase fitness in the human airway. (preprint)

SARS-CoV-2 has a broad mammalian species tropism infecting humans, cats, dogs and farmed mink. Since the start of the 2019 pandemic several reverse zoonotic outbreaks of SARS-CoV-2 have occurred in mink, one of which reinfected humans and caused a cluster of infections in Denmark. Here we investigate the molecular basis of mink and ferret adaptation and demonstrate the spike mutations Y453F, F486L, and N501T all specifically adapt SARS-CoV-2 to use mustelid ACE2. Furthermore, we risk assess these mutations and conclude mink-adapted viruses are unlikely to pose an increased threat to humans, as Y453F attenuates the virus replication in human cells and all 3 mink-adaptations have minimal antigenic impact. Finally, we show that certain SARS-CoV-2 variants emerging from circulation in humans may naturally have a greater propensity to infect mustelid hosts and therefore these species should continue to be surveyed for reverse zoonotic infections.

Mutations that adapt SARS-CoV-2 to mustelid hosts do not increase fitness in the human airway. (preprint)

SARS-CoV-2 has a broad mammalian species tropism infecting humans, cats, dogs and farmed mink. Since the start of the 2019 pandemic several reverse zoonotic outbreaks of SARS-CoV-2 have occurred in mink, one of which reinfected humans and caused a cluster of infections in Denmark. Here we investigate the molecular basis of mink and ferret adaptation and demonstrate the spike mutations Y453F, F486L, and N501T all specifically adapt SARS-CoV-2 to use mustelid ACE2. Furthermore, we risk assess these mutations and conclude mink-adapted viruses are unlikely to pose an increased threat to humans, as Y453F attenuates the virus replication in human cells and all 3 mink-adaptations have minimal antigenic impact. Finally, we show that certain SARS-CoV-2 variants emerging from circulation in humans may naturally have a greater propensity to infect mustelid hosts and therefore these species should continue to be surveyed for reverse zoonotic infections.

TMPRSS2 promotes SARS-CoV-2 evasion from NCOA7-mediated restriction (preprint)

Interferons play a critical role in regulating host immune responses to SARS-CoV-2, but the interferon (IFN)-stimulated gene (ISG) effectors that inhibit SARS-CoV-2 are not well characterized. The IFN-inducible short isoform of human nuclear receptor coactivator 7 (NCOA7) inhibits endocytic virus entry, interacts with the vacuolar ATPase, and promotes endo-lysosomal vesicle acidification and lysosomal protease activity. Here, we used ectopic expression and gene knockout to demonstrate that NCOA7 inhibits infection by SARS-CoV-2 as well as by lentivirus particles pseudotyped with SARS-CoV-2 Spike in lung epithelial cells. Infection with the highly pathogenic, SARS-CoV-1 and MERS-CoV, or seasonal, HCoV-229E and HCoV-NL63, coronavirus Spike-pseudotyped viruses was also inhibited by NCOA7. Importantly, either overexpression of TMPRSS2, which promotes plasma membrane fusion versus endosomal fusion of SARS-CoV-2, or removal of Spikes polybasic furin cleavage site rendered SARS-CoV-2 less sensitive to NCOA7 restriction. Collectively, our data indicate that furin cleavage sensitizes SARS-CoV-2 Spike to the antiviral consequences of endosomal acidification by NCOA7, and suggest that the acquisition of furin cleavage may have favoured the co-option of cell surface TMPRSS proteases as a strategy to evade the suppressive effects of IFN-induced endo-lysosomal dysregulation on virus infection.

Reduced neutralisation of the Delta (B.1.617.2) SARS-CoV-2 variant of concern following vaccination (preprint)

Vaccines are proving to be highly effective in controlling hospitalisation and deaths associated with SARS-CoV-2 infection but the emergence of viral variants with novel antigenic profiles threatens to diminish their efficacy. Assessment of the ability of sera from vaccine recipients to neutralise SARS-CoV-2 variants will inform the success of strategies for minimising COVID19 cases and the design of effective antigenic formulations. Here, we examine the sensitivity of variants of concern (VOCs) representative of the B.1.617.1 and B.1.617.2 (first associated with infections in India) and B.1.351 (first associated with infection in South Africa) lineages of SARS-CoV-2 to neutralisation by sera from individuals vaccinated with the BNT162b2 (Pfizer/BioNTech) and ChAdOx1 (Oxford/AstraZeneca) vaccines. Across all vaccinated individuals, the spike glycoproteins from B.1.617.1 and B.1.617.2 conferred reductions in neutralisation of 4.31 and 5.11-fold respectively. The reduction seen with the B.1.617.2 lineage approached that conferred by the glycoprotein from B.1.351 (South African) variant (6.29-fold reduction) that is known to be associated with reduced vaccine efficacy. Neutralising antibody titres elicited by vaccination with two doses of BNT162b2 were significantly higher than those elicited by vaccination with two doses of ChAdOx1. Fold decreases in the magnitude of neutralisation titre following two doses of BNT162b2, conferred reductions in titre of 7.77, 11.30 and 9.56-fold respectively to B.1.617.1, B.1.617.2 and B.1.351 pseudoviruses, the reduction in neutralisation of the delta variant B.1.617.2 surpassing that of B.1.1351. Fold changes in those vaccinated with two doses of ChAdOx1 were 0.69, 4.01 and 1.48 respectively. The accumulation of mutations in these VOCs, and others, demonstrate the quantifiable risk of antigenic drift and subsequent reduction in vaccine efficacy. Accordingly, booster vaccines based on updated variants are likely to be required over time to prevent productive infection. This study also suggests that two dose regimes of vaccine are required for maximal BNT162b2 and ChAdOx1-induced immunity.

A Prenylated dsRNA Sensor Protects Against Severe COVID-19 and is Absent in Horseshoe Bats (preprint)

Cell autonomous antiviral defenses can inhibit the replication of viruses and reduce transmission and disease severity. To better understand the antiviral response to SARS-CoV-2, we used interferon-stimulated gene (ISG) expression screening to reveal that OAS1, through RNase L, potently inhibits SARS-CoV-2. We show that while some people can express a prenylated OAS1 variant, that is membrane-associated and blocks SARS-CoV-2 infection, other people express a cytosolic, nonprenylated OAS1 variant which does not detect SARS-CoV-2 (determined by the splice-acceptor SNP Rs10774671). Alleles encoding nonprenylated OAS1 predominate except in people of African descent. Importantly, in hospitalized patients, expression of prenylated OAS1 was associated with protection from severe COVID-19, suggesting this antiviral defense is a major component of a protective antiviral response. Remarkably, approximately 55 million years ago, retrotransposition ablated the OAS1 prenylation signal in horseshoe bats (the presumed source of SARS-CoV-2). Thus, SARS-CoV-2 never had to adapt to evade this defense. As prenylated OAS1 is widespread in animals, the billions of people that lack a prenylated OAS1 could make humans particularly vulnerable to the spillover of coronaviruses from horseshoe bats.

In vitro evolution of Remdesivir resistance reveals genome plasticity of SARS-CoV-2 (preprint)

Remdesivir (RDV) is used widely for COVID-19 patients despite varying results in recent clinical trials. Here, we show how serially passaging SARS-CoV-2 in vitro in the presence of RDV selected for drug-resistant viral populations. We determined that the E802D mutation in the RNA-dependent RNA polymerase was sufficient to confer decreased RDV sensitivity without affecting viral fitness. Analysis of more than 200,000 sequences of globally circulating SARS-CoV-2 variants show no evidence of widespread transmission of RDV-resistant mutants. Surprisingly, we also observed changes in the Spike (i.e., H69 E484, N501, H655) corresponding to mutations identified in emerging SARS-CoV-2 variants indicating that they can arise in vitro in the absence of immune selection. This study illustrates SARS-CoV-2 genome plasticity and offers new perspectives on surveillance of viral variants.

Meta-analysis of virus-induced host gene expression reveals unique signatures of immune dysregulation induced by SARS-CoV-2 (preprint)

The clinical outcome of COVID-19 has an extreme age, genetic and comorbidity bias that is thought to be driven by an impaired immune response to SARS-CoV-2, the causative agent of the disease. The unprecedented impact of COVID-19 on global health has resulted in multiple studies generating a variety of large gene expression datasets in a relatively short period of time. In order to better understand the immune dysregulation induced by SARS-CoV-2, we carried out a meta-analysis of these transcriptomics data available in the published literature. Datasets included both those available from SARS-CoV-2 infected cell lines in vitro and those from patient samples. We focused our analysis on the identification of viral perturbed host functions as captured by co-expressed gene module analysis. Transcriptomics data from lung biopsies and nasopharyngeal samples, as opposed to those available from other clinical samples and infected cell lines, provided key signatures on the role of the host's immune response on COVID-19 pathogenesis. For example, severity of infection and patients' age are linked to the absence of stimulation of the RIG-I-like receptor signaling pathway, a known critical immediate line of defense against RNA viral infections that triggers type-I interferon responses. In addition, co-expression analysis of age-stratified transcriptional data provided evidence that signatures of key immune response pathways are perturbed in older COVID-19 patients. In particular, dysregulation of antigen-presenting components, down-regulation of cell cycle mechanisms and signatures of hyper-enriched monocytes were strongly correlated with the age of older individuals infected with SARS-CoV-2. Collectively, our meta-analysis highlights the ability of transcriptomics and gene-module analysis of aggregated datasets to aid our improved understanding of the host-specific disease mechanisms underpinning COVID-19.

Genomic epidemiology of SARS-CoV-2 spread in Scotland highlights the role of European travel in COVID-19 emergence (preprint)

SARS-CoV-2, the causative agent of COVID-19, emerged in Wuhan, China in December 2019 and spread rapidly throughout the world. Understanding the introductions of this new coronavirus in different settings may assist control efforts and the establishment of frameworks to support rapid response in future infectious disease outbreaks. We investigated the first four weeks of emergence of the SARS-CoV-2 virus in Scotland after the first case reported on the 1st March 2020. We obtained full genome sequences from 452 individuals with a laboratory-confirmed diagnosis of COVID-19, representing 20% of all cases until 1st April 2020 (n=2310). This permitted a genomic epidemiology approach to study the introductions and spread of the SARS-2 virus in Scotland. From combined phylogenetic and epidemiological analysis, we estimated at least 113 introductions of SARS-CoV-2 into Scotland during this period. Clusters containing multiple sequences suggestive of onward transmission occurred in 48/86 (56%). 42/86 (51%) clusters had no known international travel history indicating undetected introductions. The majority of viral sequences were most closely related to those circulating in other European countries, including Italy, Austria and Spain. Travel-associated introductions of SARS-CoV-2 into Scotland predated travel restrictions in the UK and other European countries. The first local transmission occurred three days after the first case. A shift from travel-associated to sustained community transmission was apparent after only 11 days. Undetected introductions occurred prior to the first known case of COVID-19. Earlier travel restrictions and quarantine measures might have resulted in fewer introductions into Scotland, thereby reducing the number of cases and the subsequent burden on health services. The high number of introductions and transmission rates were likely to have impacted on national contact tracing efforts. Our results also demonstrate that local real-time genomic epidemiology can be used to monitor transmission clusters and facilitate control efforts to restrict the spread of COVID-19. FundingMRC (MC UU 1201412), UKRI/Wellcome (COG-UK), Wellcome Trust Collaborator Award (206298/Z/17/Z - ARTIC Network; TCW Wellcome Trust Award 204802/Z/16/Z Research in contextO_ST_ABSEvidence before this studyC_ST_ABSCoronavirus disease-2019 (COVID-19) was first diagnosed in Scotland on the 1st of March 2020 following the emergence of the causative severe acute respiratory system coronavirus 2 (SARS-CoV-2) virus in China in December 2019. During the first month of the outbreak in Scotland, 2310 positive cases of COVID-19 were detected, associated with 1832 hospital admissions, 207 intensive care admissions and 126 deaths. The number of introductions into Scotland and the source of those introductions was not known prior to this study. Added value of this studyUsing a combined phylogenetic and epidemiological approach following real-time next generation sequencing of 452 SARS-CoV-2 samples, it was estimated that the virus was introduced to Scotland on at least 113 occasions, mostly from other European countries, including Italy, Austria and Spain. Localised outbreaks occurred in the community across multiple Scottish health boards, within healthcare facilities and an international conference and community transmission was established rapidly, before local and international lockdown measures were introduced.